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Abstract 

We focus on the transition from quantum to classical behavior in thermodynamic functions 
and time correlation functions of a system consisting of three identical quantum spins s 
that interact via isotropic Heisenberg exchange. The partition function and the zero-field 
magnetic susceptibility are readily shown to adopt their classical forms with increasing s. 
The behavior of the spin autocorrelation function (ACF) is more subtle. Unlike the classical 
Heisenberg trimer where the ACF approaches a unique non-zero limit for long times, for 
the quantum trimer the ACF is periodic in time. We present exact values of the time average 
over one period of the quantum trimer for s < 7 and for infinite temperature. These aver- 
ages differ from the long-time limit, (9/40) In 3 + (7/30), of the corresponding classical 
trimer by terms of order 1/s^. However, upon applying the Levin n-sequence accelera- 
tion method to our quantum results we can reproduce the classical value to six significant 
figures. 
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1 Introduction and summary 



In recent years considerable attention has been devoted to the magnetic properties 
of synthesized organic complexes ("molecular magnets") containing small numbers 
of paramagnetic ions [1-5]. With the ability to control the placement of magnetic 
moments of diverse species within stable molecular structures, it is now possible to 
test basic questions concerning their magnetic properties and to explore the design 
of novel systems that offer the prospect of useful applications [6,7]. Intermolecular 
magnetic interactions are typically extremely weak compared to intramolecular in- 
teractions, so a bulk sample can be considered as independent individual molecular 
magnets. 

The present study is motivated by the successful synthesis of two trimers, one 
[8] consisting of V^^ ions (s = 1/2) and the second [9] consisting of Fe^+ ions 
(s = 5/2). With the anticipation of the successful synthesis of yet other trimers, in 
this article we consider the thermodynamic functions and the time correlation func- 
tions for three identical quantum spins s which interact via isotropic Heisenberg 
exchange in the absence of an external magnetic field. We are especially interested 
in comparing the behavior of these quantities for diverse s, and in particular in the 
transition to classical behavior which occurs for large s. This transition is readily 
analyzed for the thermodynamic functions. In particular, we show that the partition 
function and the zero-field magnetic susceptibility adopt the correct classical forms 
[10] with increasing s. Our analysis also provides results for the deviation from 
classical behavior depending on the size of s and the temperature of the system. 

Knowledge of the exact two-spin time correlation functions is of great importance 
since this enables one to predict the outcome of measurements such as the proton- 
spin lattice relaxation rate [11] by nuclear magnetic resonance (NMR) techniques 
and inelastic magnetic neutron scattering [12]. In an earlier publication [13] we ex- 
amined in detail the properties of the thermal equilibrium autocorrelation function, 
to be denoted by As{t, T), for a quantum Heisenberg dimer composed of two iden- 
tical spins s which interact with isotropic exchange. Some features of As{t, T) that 
occur for a dimer are readily shown to occur for the trimer, too, and we will not 
deal with those issues. However, there are two features of ^^(t, T) for the quantum 
trimer which deserve careful attention and these are considered in detail. 

First, for the quantum Heisenberg trimers the quantity As{t, T) is a periodic func- 
tion of the time with a recurrence time r. For the corresponding classical Heisen- 
berg trimer, whose ACF is denoted by Ac{t,T), there exists a unique, non-zero 
long-time limit Ac{oo,T) [14]. To make a meaningful comparison between the 
asymptotic, long-time behavior of the classical Heisenberg trimer with a quantum 
trimer, we suggest comparing Ac(oo, T) with the time average of As{t, T) over the 
corresponding recurrence time r, to be denoted by Ag (T) ; this value is equivalent to 
the coefficient of 5{uj) in the expression for the Fourier time transform of As{t, T). 



In particular we will explore the large s behavior of As{T). Second, for the classical 
Heisenberg trimer, and for the special case of infinite temperature it has been shown 
that Ac{oo, oo) = (9/40) In 3 + 7/30 [14]. This curious exact value differs from the 
result, 1/N, that follows from a phenomenological diffusive spin dynamics (DSD) 
[14] based on linear equations of motion for a ring of N, here A^ = 3, interacting 
classical Heisenberg spins. By contrast, for both the quantum and classical dimers 
it has been found that As{oo) = Ac{oo,oo) = |, independent of s and in agree- 
ment with the DSD result. It is thus of interest to track the emergence of the exact 
classical result upon considering the sequence of quantum trimers for increasing s. 

We have derived the exact values of ^^(oo) for the quantum trimers for the choices 
s = 1/2, 1, 3/2, . . . , 7. It is necessary to consider the results for half-integer s 
separately from those derived for integer s since they exhibit different behaviour 
patterns. Each subsequence appears to converge very slowly to the above value of 
Ac{oo,oo) for the classical trimer. For spins s the deviation of ^^(oo) from the 
classical result is of order 1/s^. However, when we apply the Levin w-sequence 
acceleration method [15,16] to the subsequence for half-integer s we arrive at an 
estimate for the large-s-limit which agrees to 6 significant figures with Ac{oo, oo). 
For integer s the Levin w-estimate agrees with Ac(oo, oo) to 5 significant figures. 

The layout of this paper is as follows. In Sec. 2 we calculate the partition function 
and the zero-field susceptibility of the quantum trimer for general s. We show that 
both quantities approach the corresponding results for the classical trimer. In Sec. 3 
we derive a general formula for the time average As{oo). Numerical evaluation of 
the formula becomes a very lengthy process with increasing s, however, we have 
been able to perform these calculations for s < 7. The Levin w-estimates for A^ (oo) 
are also provided in Sec. 3. 



2 Thermodynamic functions 

The quantum trimer is specified by the Hamilton operator 

^°^n2 fti ■ I2 + I2 ■ Is + Is ■ |i) (1) 

= ^ i§ - £1 - £2 - £s j ; § = Sl + S2 + £3 , 

where the spin operators satisfy the usual commutation relations and where J has 
units of energy. J > describes antiferromagnetic and J < ferromagnetic cou- 
pling. Throughout this article it is assumed that the spin quantum numbers of the 
three sites are identical, Si = S2 = s-^ = s. The eigenstates | S", M, 6*23 ) of the 
Hamilton operator can be chosen as simultaneous eigenstates of the total spin S"^, 



its z-component Sz, and of 523 ~ (I2 + Is)^- The quantum numbers S, M, S23 
characterize the eigenstates completely. The eigenvalues Es of the Hamilton oper- 
ator 



Es = ^ iSiS + l)-3sis + l)) (2) 

do not depend on S23 or M. Thus the partition function in the canonical ensemble 
reads 



S,M,S23 

The respective classical Hamilton function He is defined as [13] 



Hc = Jc (ei-e2 + 62 -63 + 63-61) , Jc = J s{s + I) , (4) 

where ei, 62 and 63 are unit vectors (c-numbers). Then the classical partition func- 
tion turns out to be [10] 



Z, = I dED,iE)exp{-PE), (5) 

with the classical density of states 



DciE) 



2fc\/f + 3 


-lJc<E< -J, 




V 


-Jc< E < 3Jc 
else 



(6) 



which is normalized to unity. In order to compare quantum and classical density of 
states, the energy spectra of the quantum trimers for different s have to be mapped 
onto the same energy interval; we take [— 3Jc/2, 3 JJ, i.e. all energies are divided 

bys(s + l). 

Figure 1 demonstrates how the quantum density of states approaches the classical 
limit with increasing spin s. Considering the eigenvalues (2) of the Hamilton opera- 
tor (1) and their multiplicities, which originate from the degeneracy in the magnetic 
quantum number M and from different possibilities to couple to a certain total spin 
S, one can show analytically that the quantum density of states converges against 
the classical one for infinite s [17]. 

But the coincidence of the classical and the quantum density of states for high spin 
does not mean that other observables coincide, too. An interesting example is given 
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Fig. 1. Normalized density of states for quantum Heisenberg trimers (dashed lines) and 
their classical counterpart (solid line). 

by the zero field susceptibility. The susceptibility is defined as the derivative of the 
magnetisation 



M = ^tT {-gfiBS.0 ^-} , H = Ho + g^sBS, 



with respect to the magnetic field B 



(7) 



Xo^ 



'dM 



9 f^sP 



B=0 



trOS 



-I3H 



B=0 



(8) 



For ferromagnetic coupling, Fig. 2 r.h.s., the graphs for different s nearly coin- 
cide with each other and with the classical result. However, for antiferromagnetic 
coupling the zero field susceptibility behaves differently for integer and half- integer 
spin quantum numbers as shown on the l.h.s. of Fig. 2. The reason is that for integer 
spin the ground state, which is non-degenerate, has S = and thus the zero field 
susceptibility approaches zero for small temperatures T, whereas for half-integer 
spins the ground state, which is fourfold degenerate, has S = 1/2 which causes the 
susceptibility to go to infinity for small T. However, if we consider a fixed nonzero 
value of T, the susceptibilities tend to the classical result with increasing s. The 
classical result is "indifferent" at T = 0, it approaches 1. 



3 Autocorrelation function 



Another important observable is the two-spin time correlation function because it 
serves as a major ingredient for several quantities such as the spin lattice relaxation 
rate and the neutron scattering cross section. 
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Fig. 2. Zero-field susceptibility for A^ = 3 and choices s = 1/2, 1, . . . , 9. The solid lines 
show the resuh for half integer spins, the dashed lines for integer spin quantum numbers. 
The classical result is given by the thick dotted line. L.h.s.: antiferromagnetic coupling, 
r.h.s.: ferromagnetic coupling. 

Considering that the Hamilton operator (1) is isotropic, one obtains for the autocor- 
relation function 



^^ Re {Es,m,sJS, M, S2S I si.(t) ■ s,M e'^^" I S, M, S^s )] 



/3Ho 



Es,M,sJS, M, ^23 I £1.(0) ■ £1.(0) e-"^' I S, M, S23 ) 



Since ( S, M, ^23 I £1. I S', M', S'^^ ) is zero if |5 - 5'| > 1, the contributing fre- 
quencies are a; = J{S + l)/h. They are all multiples of a basic frequency, which is 
cUr = J/h for integer spin quantum numbers and cUr = J/{2h) for half integer spin 
quantum numbers. Thus the autocorrelation function is periodic with a recurrence 
time 2tx jijjr- 

We have evaluated the expression ^^(oo) 



T~*^Jq 



As{oo)= lim 1" dt As{t,T) 



(10) 



for half-integer values 1/2 < s < 13/2 and for integer values 1 < s < 7, and these 
are listed in table 1. It appears to be impractical to extend these results to larger 
spin values since the amount of computer time required grows at an astonishing rate 
with increasing s. Fortunately additional calculations are unwarranted. In Fig. 3 we 
display our results versus the independent variable l/(2s + 1)^ along with the solid 
line which has been chosen to pass through the quantum results for s=ll/2 and 
s=13/2. The good agreement between the results for the larger half- integer values 
of s and the solid line is consistent with the conclusion that the deviation between 
the quantum results and the classical trimer decreases monotonically to zero but 
very slowly, the deviation being of order l/(2s + 1)^. In fact, if we approximate 
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Fig. 3. Results for ^^(oo) (open and full circles) and Ac(oo, oo) (cross). The line, which 
connects the results for s = 11/2 and s = 13/2, is drawn to guide the eye. 

^^(oo) by the form 



AJoo) 



A 



B 



(2s + 1) 



(11) 



we may use our results for s = 11/2 and 13/2 to determine the unknown param- 
eters A and B. In particular the result A = 0.480511 provides an estimate for 
lims-,QoAs{oo). This result is rather close to the exact classical result [14] 



Ae(oo, oo) = (9/40) ln3 + 7/30 = 0.4805210983 . 



(12) 



Adopting Eq. (1 1) for the case of integer values of s and using our results for s = 6 
and s = 7 we find that A = 0.480575. 

We may obtain an improved estimate for lims^QoAs{oo) by exploiting the Levin 
M-sequence acceleration method [15,16] which is tailor-made for such slowly con- 
vergent, monotonic sequences as the ones we face. If the sequence elements are 



labelled Ui, U2, U3, . . . (to be identified with ^41/2(00), ^3/2(00), . . . ) and if we 
define the quantities ui = Ui,U2 = U2 — Ui,u^ = U3 — U2, ■ ■ ■ then the Levin 
M-estimate for lirrin^ooUn based on employing the first M values of f/„ is given by 



U[M] 






(13) 



We find that U[7] = 0.48052085 .... For the corresponding sequence of integer 
values of s we find that U[7] = 0.4805179 .... 

It is interesting to note that for half-integer values of s the Levin w-estimates are 
closer to the exact classical result than those for integer s. 
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.49382716 

0.48648889 

0.48385973 

0.48264642 

0.48199074 

= 0.48159723 

= 0.48134282 



Table 1 

Time average, ^^(oo), of the quantum autocorrelation function in the high-temperature 

limit for s = 1/2, 1,..., 7. 
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